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We calculate the spin temperature and 21 cm brightness of early Hll regions 
around the first stars. We use outputs from cosmological radiation-hydrodynamics 
simulations of the formation and evolution of early Hll regions. In the pre-reionization 
era, Hn regions around massive primordial stars have diameters of a few kpc. The gas 
within the Hll regions is almost fully ionized, but begins recombining after the central 
stars die off. The relic Hn regions are then seen as bright emission sources in hydrogen 
21 cm. We make brightness temperature maps of the Hn regions, accounting for radia- 
tive coupling with Lyman-a photons in a simplified manner. The spin temperature in 
the relic Hn region is close to the gas kinetic temperature, generally several hundred 
to several thousand degrees. We show that the relic Hn region can be as bright as 
6Tb ~ 100 mK in differential temperature against the cosmic microwave background 
for an angular resolution of sub-arcseconds. While individual early Hn patches will 
not be identified by currently planned radio telescopes, the collective fluctuations from 
early Hn regions might imprint signatures in the 21 cm background. 
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1 INTRODUCTION 

Reionization of hydrogen in the intergalactic medium (IGM) 
is an important landmark in the history of the universe. It 
is generally thought that the sources of ionizing photons are 
early generations of stars and galaxies formed at z > 6. 
Recent observations provided important constraints on the 
epoch of reionization; the Gunn-Peterson troughs are found 
in the spectra of distant quasars at z = 6 — 7 (White et 
al. 2003; Fan et al. 2006), while the large-scale polarization 
of cosmic microwave background (CMB) measured by the 
WMAP satellite suggests that reionization may have begun 
at z > 10 (Page et al. 2007; Komatsu et al. 2008). 

Observing 21 cm emission or absorption from hydrogen 
atoms in the high redshift IGM is a promising way of reveal- 
ing the detailed process of reionization. There have been a 
number of studies on 21 cm signature of large-scale reion- 
ization by galaxies (McQuinn et al. 2006; Iliev et al. 2006; 
see Furlanetto et al. 2006 for a review). These studies are 
aimed at making forecasts for near-future radio observations 
such as LOFAR, and thus consider large-scale structures of 
ionized/neutral IGM at z ~ 6 — 10. 

There have been also theoretical studies on 21 cm sig- 
natures from pre-reionization epochs. Kuhlen et al (2006) 
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study the 21 cm signature of early mini-quasars. They show 
that X-rays from the mini-quasar raise the gas kinetic tem- 
perature and enhance 21 cm signals. Chen & Miralda-Escude 
(2004; 2008) argue that X-rays from the first stars heat the 
surrounding gas and couple the spin temperature to its ki- 
netic temperature, generating a large Lyman-a absorption 
sphere. Shapiro et al (2006) estimate the 21 cm fluctuations 
caused by cosmological "minihalos" and by the IGM. They 
conclude that the 21 cm emission from minihalos dominates 
over that from the diffuse IGM at z < 20. Furlanetto & Oh 
(2006) argue, however, that minihalos generate only small 
21 cm fluctuations. Clearly, it is important to identify dom- 
inant sources of 21 cm fluctuations at high-redshifts. None 
of these works, however, consider 21 cm signals from early 
relic Hll regions, which are hot and (partially) ionized even 
after the central sources' lifetimes, and hence can be bigger 
and brighter than cosmological mini-halos. Interestingly, in 
the standard theory of cosmic structure formation based on 
cold dark matter, there is a large gap between the time when 
the first stars are formed (Tegmark et al. 1997; Yoshida et 
al. 2003; 2008) and when the IGM is completely ionized. 
Hence it is expected that there were early relic Hll regions 
during the reionization process. Although individual Hll re- 
gions are too small to be detected even by future radio tele- 
scopes, their collective signals may imprint distinguishable 
fluctuations. In this paper, we study the 21 cm signature of 
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early relic Hn regions around the first stars using cosmolog- 
ical simulations. We assume that massive primordial stars 
are formed in early low-mass dark matter halos and ionize a 
large volume of the surrounding IGM in their short lifetime 
of a few million years. We compute the spin temperature 
and 21 cm differential temperature fluctuations of the relic 
Hn regions. 



2 NUMERICAL SIMULATIONS 

We use the outputs of a cosmological simulation of Yoshida 
et al. (2007, hereafter Y07). Briefly, they carried out three- 
dimensional radiation hydrodynamics calculations of the for- 
mation and evolution of early Hn / Hem regions around the 
first stars. Primordial gas chemistry such as hydrogen re- 
combination is followed in a non-equilibrium manner in the 
simulation. We generate the same realization as that of Y07 
but for a slightly lower value of the normalization of the 
density fluctuation amplitude, ag = 0.8. 

The other cosmological parameters are the same as 
in Y07, (On, Ov, Ho) = (0.26, 0.04, 0.7, 70[km- S - 1 • 
Mpc -1 ]), which are consistent with the standard ACDM 
model (Komatsu et al. 2008). Lowering the value of ag sim- 
ply shifts the formation epoch of the first star to a later 
time. The main properties of the Hn region are almost the 
same as those presented in Y07. In our simulation, the first 
star-forming cloud is located at z = 22. We embed a Popu- 
lation III star with a mass of WOMq. In order to focus on 
the evolution of an early Hn region, we assume that the cen- 
tral star dies without triggering supernova explosion. Such 
a massive star may also leave behind a remnant black hole. 
We defer detailed studies on the feedback effects from su- 
pernovae and from early remnant blackholes to future work. 
We note that early Population III remnant blackholes do 
not accrete the surrounding gas and emit X-rays efficiently, 
due to photoevaporation of gas from the host halo by the 
progenitor star, as well as photo-heating by the accreting 
black hole itself (Y07; Johnson & Bromm 2007; Alvarez et 
al. 2008). We further assume that UV radiation background 
and UV radiation from other sources do not exist and that 
the central star is the first radiation source. We follow the 
thermal evolution of the relic Hn region for over one hundred 
million years under these assumptions. 



3 21 CM SIGNATURES FROM EARLY Hn 
REGIONS 

The 21 cm line is produced via the transition between the 
triplet and the singlet sublevels of the hyperfine structure 
of the ground level of a neutral hydrogen atom. This wave- 
length corresponds to a frequency of 1420 MHz and a tem- 
perature of T* = 0.068K. The transition between the levels 
occurs through three processes: absorption of CMB photons, 
collisions with other particles (hydrogen atoms, free elec- 
trons and protons), and scattering of UV photons. 

The spin temperature is expressed as a weighted mean 
of these processes as (Field 1958) 



where Tk is the gas kinetic temperature, T c is the color tem- 
perature of Lyman-a photons, and x c and coefficients 
for collisional and radiative coupling, respectively. 

We do not consider a cosmological UV background ra- 
diation because the 21 cm differential brightness tempera- 
ture of relic Hn regions is not significantly affected by the 
strength of radiative coupling. Also, we do not include the ef- 
fect of heating by X-ray photons from the first stars. The gas 
outside the Hn region can be significantly heated by X-ray 
induced secondary electrons. However, for a 100 Mq star, 
the size of the heated region is only slightly larger than the 
extent of the Hn region (Chen & Miralda-Escude 2008, Fig- 
ures 3-5). Hence ignoring the X-ray heating effect does not 
affect significantly our results; in this case it would boost 
the observable signal by a factor of order unity. For very 
massive stars (> 300Mq), which are hotter and generate 
more X-rays, the signal from the X-ray heated region could 
be significantly larger. 

The coupling coefficient for collisions, x c , is given by 



A 10 T K 



(C h + C c + C p ) 



(2) 



where the spontaneous emission rate Aio is 2.85 x 10 _15 s _1 , 
and Ch, C e and C p are the de-excitation rates for collisions 
with neutral atoms, free electrons and protons, respectively 
0. We use fitting functions of Zygelman (2005) and Liszt 
(2001) for H-H collisions and e-H collisions, respectively. For 
p-H collisions, we use the coefficient 3.2 times as large as that 
for H-H collisions (Smith 1966). 

The radiative coupling included in equation JTJ is 
known as the Wouthuysen-Field effect. It is extremely dif- 
ficult to evaluate accurately the coupling coefficient for the 
incident UV radiation, x a , because one needs to solve, in 
principle, full radiation transfer equations including the ef- 
fect of absorption and re-emission of Lyman-a photons, the 
Hubble flow, and the injection of new photons. By noting 
that tight radiative coupling is easily achieved under many 
cosmological conditions (e.g., Furlanetto et al. 2006), we as- 
sume that Lyman-a photons produced by recombination is 
sufficient for effective radiative coupling, and examine the 
overall effect by setting a constant x a = 0.5 within Hn re- 
gions. We compute the color temperature, T c , by solving the 
following equation 



1 + T 5 c/Tk 
1 + T sc /T s 



(3) 



Is = 



■ x c Tk + x a T c 



1 -(~ Xc ~t~ Xa 



(i) 



where T sc = (2/9)Tk^fi/AfD,a> which is obtained by 
solving the radiative transfer equation for scattering of 
UV photons for large optical depths (Furlanetto et al. 
2006). Here, V21 is the frequency of 21 cm line, 
^/(2fc B T K )/(mHC 2 ) v a is the Doppler width of Lyman-a, and 
v a is the central Lyman-a frequency (2.47 x 10 15 Hz). We 
determine the color temperature and the spin temperature 
by solving equations ([TJ and <(3j iteratively. 

We briefly describe the evolution of the Hn region. 
Shortly after the central star dies, the kinetic temperature 
in the Hn region is higher than several thousand Kelvin. 



Although the simulation of Y07 includes He atoms and ions, 
we ignore the collisions between H atoms and He species because 
the number fraction of He is small. 
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Figure 1. Wc plot the 21 cm spin temperature against gas den- 
sity. Left: No radiative coupling. Right:With radiative coupling 
for x a = 0.5. We use an output at z = 21.7 when Tqmb = 62 K. 

While high density regions within the Hll region recombine 
and cool quickly, a large fraction of the Hll region has a low 
density (n < 0.01cm -3 ), and thus recombination and gas 
cooling occur rather slowly. The gas temperature remains a 
few thousand Kelvin at 30 Myrs, and the ionization fraction 
is down to about 10%, 5% at 30 Myrs and 50 Myrs after the 
star dies, respectively. 

Figure Q] shows the 21 cm spin temperature against gas 
density at a time shortly after the ionizing source has been 
turned off. To see the effect of radiative coupling, we calcu- 
late the spin temperature ignoring (left panel) and account- 
ing for (right panel) radiative coupling. At low densities, 
the spin temperatures are close to the CMB temperature 
for both cases, because the collisional coupling is weak. The 
gas within the Hll regions has been once photo-ionized, and 
so the gas kinetic temperature is high, which brings the spin 
temperature above ~ 100 K for the case without radiative 
coupling (left) and above ~ 10 3 K with radiative coupling 
(right). Clearly, radiative coupling raises the spin tempera- 
tures of intermediate-density, ionized gases. We see another 
branch in the high density, low spin temperature portion in 
this phase diagram. It corresponds to the gas outside the 
relic Hll region, which has been neutral and hence has a low 
kinetic temperature. Note that there are also neutral regions 
where the spin temperatures are slightly lower than CMB 
temperature because of effective collisional coupling. 

We now observe the simulated Hll region in 21 cm. The 
radiative transfer equation for 21 cm radiation for brightness 
temperature, Tb, is given by 

T b = T C MBe- T + T s (l-e- T ). (4) 

In practice, we discretize the equation as in Kuhlen et al. 
(2006) to calculate the brightness temperature for our sim- 
ulated Hll region. The observed flux of 21 cm line can be 
expressed by the differential brightness temperature against 
the CMB temperature as ST b = (T b - Tcmb)/(1 + z). To 
make maps of Ts and 5Tb, we use a cubic region of 100/i _1 
comoving kpc on a side, centered at the Hn region, and then 
divide the box into 128 3 cells. Then the angular resolution 
of the map is about 0.02 arcseconds. We use the central 
slice of with a 7.8/i -1 comoving kpc depth along the line 
of sight to make the maps shown in Figure [2] The corre- 
sponding spectral width is 190 Hz at z=20. Note that 




Figure 2. We show the evolution of the spin temperature. The 
output time of the top panels is z = 21.7, when the ionizing source 
is switched off, and that of the bottom panels is z = 19.5, about 
30 Myrs later. Left: No radiative coupling. Right: With radiative 
coupling (x a = 0.5). The plotted region is 100 h^ 1 comoving kpc 
on a side. 



the required angular and spectral resolution are out of reach 
of any planned instruments; however, the collective effect of 
such relic Hn regions might be detectable. 

Figure [2] shows the mass-weighted mean spin temper- 
ature at two output times, when the ionizing source died 
(top panels) and about 30 Myrs after (bottom panels). The 
bright, colored area is almost fully ionized (i.e. Hn region) 
and the radius of the Hll region is about 50/i -1 comoving 
kpc. At the top panels, the spin temperature ranges from a 
few tens Kelvin (~ Tomb) to several thousand Kelvin. The 
top-right panel in Figure[2]shows that the radiative coupling 
significantly affects the spin temperature in the Hll region. 
After the central star died (t — 0), the Hll region starts re- 
combining, and the gas kinetic temperature decreases. Over- 
all, the mean spin temperature at t = 30 Myrs is smaller 
than at t = 0. Hydrodynamic effects are also seen in this 
high-resolution map. Strong Shockwaves were driven during 
the early evolution of the Hll region (Kitayama et al. 2004; 
Y07), which are seen as a bright "ring" in the bottom panels 
in Figured] 

Figure shows the time evolution of the resulting 21 
cm signal for the case with no radiative coupling. Regions 
where ST^ are colored in black. While a large volume 
surrounding the Hll region has T sp m ~ Tcmb (see Figure [2}, 
which appear in absorption. We also calculated the differen- 
tial temperature for the case with efficient radiative coupling 
and find little difference with Figure [3] The clear difference 
in the spin temperature shown in Figure [2] does not affect 
the differential temperature. This is because, as long as Ts 
is much higher than Tcmb , the observed brightness is deter- 
mined by the column density of neutral hydrogen and the 
width of 21 cm line, which are independent of Ts (e.g., Scott 
& Rees 1990). The large filamentary structure seen in Figure 
[3]is a neutral region which was not completely ionized by the 
central star because of its strong self-shielding. On average, 
its spin temperature is not high, hence the structure is not 
prominent in Figure [2] but it shows up clearly in differential 
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Figure 3. We plot the time evolution of the 21 cm differential 
brightness temperature. Top-left:Signals at the time of the ion- 
izing source dead (t = 0; z = 21.7). Top-right, bottom-left, and 
bottom-right are signals at t = 10 Myrs, 30 Myrs, and 50 Myrs. 



temperature because it has a large column density of neutral 
hydrogen (and hence a large 21 cm optical depth). 

Figure [3] shows that the relic Hn region is a 21 cm emis- 
sion source. In the relic Hn region, STb is as high as several 
hundred mK, while in the surrounding neutral regions <5Tb 
can be as low as ~ — 10 mK. The relic Hn region is bright for 
over 50 Myrs after the source died, because the number den- 
sity of neutral hydrogen increases by recombination in the 
relic Hn region, while the spin temperature remains higher 
than Tomb- The area-averaged 8T^ remains a few tens mK 
for over 100 Myrs (at the last output of our simulations). 



4 DISCUSSION 

We have studied the 21 cm signature from relic Hn re- 
gions at high redshifts. We have computed the spin tem- 
perature and 21 cm brightness temperature using radiation- 
hydrodynamics simulations. We have examined the effect 
of radiative coupling on the spin temperature under simpli- 
fied assumptions. We have shown that, whereas the radiative 
coupling significantly affects the spin temperature, it has lit- 
tle effect on the 21 cm differential temperature. The relic Hn 
region is seen as a bright emission source in 21 cm, with the 
differential temperature being up to ~ 100 mK for very high 
angular resolutions. Previous works estimated the spin tem- 
perature and 21 cm signals of Hn regions, but they did not 
follow evolution by computing recombination, cooling and 
the expansion of Hn regions dynamically. We compared our 
results with those of Nusser (2005), who calculated the 21 
cm signatures under the isobaric condition. The overall evo- 
lution is similar, but the Hn region in our simulation shine 
in 21 cm for a longer time because recombination occurs 
slowly in low density regions. 

Early Hn regions are generally very small and thus the 
individual 21 cm sources will not be detected by the Low Fre- 
quency Array (LOFARjf) and the Mileura Widefield Array 
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(mwaJE even by the next-generation low-frequency arrays 
such as the Square Kilometre Array (SKA j3- Minihalos simi- 
larly contain little mass, and are not individually detectable. 
A very large, percolated Hn regions or a group of nearby Hn 
regions may be detectable through the effect of strong grav- 
itational lensing (Li et al. 2007). As we have shown, relic 
Hn regions are emission sources for a long time (> 50 Myrs, 
see Figure [2}, and thus the number density of such 21 cm 
sources at a given frequency can be larger than Lyman-a 
spheres around short-lived stars. Therefore, early relic Hn 
regions could imprint distinct, strong collective signatures 
in the 21 cm background. In our future work, we will calcu- 
late the abundance and the clustering of the relic Hn regions 
and derive the 21 cm fluctuation amplitudes. 

Finally, we comment that there are some other sources 
of 21 cm emission and absorption at high redshifts. Ripa- 
monti et al (2008) study the effect of X-rays from early black 
holes. They show that the differential brightness tempera- 
ture is ~ 20 - 30 mK at z ^ 12. Zaroubi et al (2007) show 
that X-ray heating produces a differential brightness tem- 
perature of the order of ~ 5 — 10 mK out to a few comoving 
Mpc distance from black holes. These brightness temper- 
atures are smaller, at least locally, than that of relic Hn 
regions studied in the present paper. Thomas & Zaroubi 
(2008) consider both primordial black holes and Population 
III stars. They argue that the heating patterns around these 
objects is significantly different. A few other exotic models 
of ionization such as ultra-high energy cosmic rays and de- 
caying dark matter particles are proposed (Shchekinov & 
Vasiliev 2007). Such different sources can likely be distin- 
guished by their different fluctuation spectra. 
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